Generalized diffusion equation 
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Modern analyses of diffusion processes have proposed nonlinear versions of the Fokker-Planck 
equation to account for non-classical diffusion. These nonlinear equations are usually constructed on 
a phenomenological basis. Here we introduce a nonlinear transformation by defining the g-generating 
function which, when applied to the intermediate scattering function of classical statistical mechan- 
ics, yields, in a mathematically systematic derivation, a generalized form of the advection-diffusion 
equation in Fourier space. Its solutions are discussed and suggest that the g-generating function 
approach should be a useful tool to generalize classical diffusive transport formulations. 
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The development of generalizations of classical statis- 
tical mechanical paradigms to describe non-equilibrium 
systems has become an important area of activity. One 
of the best known examples of this approach is the so- 
called " nonextensive" thermodynamics first introduced 
by Tsallis which is based on a new definition of entropy 
that as being the g-logarithm of the number of states, 
rather than the natural logarithm in Boltzmann's defi- 
nition. It has been shown that thermodynamics can be 
consistantly built using the Tsallis entropy as a starting 
point 0- The success of this program has led to the in- 
troduction of other "nonextensive " generalizations such 
as the nonextensive diffusion equation Q . As in the case 
of the entropy, these generalizations are usually based 
on the replacement of a classical, extensive (frequently 
linear) phenomenological expression by a parameterized, 
nonextensive (frequently non-linear) expression which re- 
duces to the classical expression for some value of the pa- 
rameter. In this paper, we show that for the particular 
case of diffusion it is possible to introduce these general- 
izations at a more fundamental level. The key point is 
that statistical mechanical derivations often make use of 
resummations and ansatzs which can be made "nonex- 
tensive " in a natural way since they are intrinsically 
heuristic. 

Processes related to diffusive transport are most com- 
monly measured by scattering methods (photon corre- 
lation and thermal neutron spectroscopy) and the spec- 
tra are analyzed using the intermediate scattering func- 
tion II 



i^,(fc,i) = (exp[zk.(R(t)-R(0))]), 



(1) 



which is the space-Fourier transform of the single parti- 
cle density correlation function, i.e. the probability that 
a tagged particle be at position R(t) at time t, given it 
was at R(0) at the initial time. In classical statistical 
mechanics, one shows that this function obeys the usual 
diffusion - or advection-diffusion - equation in the hydro- 
dynamic limit, i.e. for long wavelengths: k£ << 1, where 
£ is the mean free path. This is easily obtained, for in- 



stance, by the cumulant expansion of Fs{k,t) 4j. Here 
we generalize the procedure by defining the q-generating 
function 



(ikp 



K,(i) 



(2) 



or equivalently the g-cumulant expansion 

{thy 



Fsik,t) 



exp 



1 + (1-9)E 



(3) 



where Kj(t) are time dependent functions which can be 
expressed in terms of the spatial moments . Note that 
in the limit g — )■ 1, the q-log and the g-exp reduce to the 
usual logarithmic and exponential expressions H. From 
0, we have 



dF,ik,t) 
dt 



j=o 



dt 



FHKt), (4) 



which is an exact equation that follows from the nonlinear 
transformation jSjl. Low k expansion on the r.h.s. of 
yields 



dFs{k,t) 
dt 

where 
dKi{t) 



dKi{t) lj^^ dK2{t) 

dt 2 dt 



(5) 



- -([R(i)-R(0)]) = u+(v(i)) = u, (6) 



with u the macroscopic velocity or " drift term " , and 



dK2it) d 



dt 



dt 



([R(t)-R(0)n =2D, 



(7) 



where, for simplicity, we have considered isotropic diffu- 
sion, with D the self-diffusion coefficient. So to second 



2 



order, as usually considered for the classical hydrody- 
namic limit, we obtain the nonlinear advection-difFusion 
equation 



dFs{k,t) 
dt 



ik-uF^{k,t) - Dk^ F^{k,t), (8) 



or 



1 dF}-i{k,t) 
1-q dt 



= tk u- Dk^ 



(9) 



The solution for initial condition Fs{k,t = 0) = 1 (i.e. a 
real space delta-function) and for positive times is 

F,{k,t) = [l + {I- q){ik-u- Dk^)t]^ 

= expg [{ik ■ u - Dk^)t] , (10) 

a form commonly found in nonextensive statistical me- 
chanics 1^. It is clear that, in the limit q 1, one 
retrieves the classical result: Fs{k, t) 
Equation (O rewritten as 



^(jk-u-_Dfc'^)t 



dFs{k,t) 
dt 



-ik-Ug{k,t)Fs{k,t) = - Dk"^ F^{k,t) , (11) 



with Uq{k,t) ~ uF^^^{k,t), can be viewed as the ana- 
logue in Fourier space of the porous media equation 0,0, 
but in contrast to the latter, the nonlinear advection- 
diffusion equation ffTT^ yields generalized, but classical 
diffusive behavior in real space (i.e. with r/V Dt scaling, 
as we shall see below) , and therefore applies to a different 
class of diffusive transport. 

In the absence of flow (u = 0), we call Eq.JHJ the q- 
diffusion equation 



dF,{k,t) 
dt 



= -Dk^F^{k,t). 



(12) 



This generalized equation can be cast into a regular diffu- 
sion equation with an effective diffusion coefficient -De// 
by using the identification 



D 



D, - DFr\k,t), 



(13) 



or, substituting the solution given in IjlUI) (with u = 0) 
and expanding in Dk^t (i.e. for long wavelengths). 



(14) 



that is Dq = D + higher moments {q^ k,t), where the 
higher moments depend on q, and also on k and t as for 
the space- and time-dependent transport coefficients in 
generalized hydrodynamics fij. For q — 1, only the term 
J = survives in the sum in (|14|l . and one retrieves the 
classical diffusion formulation. 

We now consider what type of distributions are ob- 
tained from the solution to the g-diffusion equation (fT^ : 



Fs{k,t) = [l-{l-q)Dkh] 



(15) 



As the scattering function must be real and non-negative, 
the argument in brackets must be positive. This is always 
realized for g > 1 . If we are to allow g < 1 , then it is nec- 
essary to limit the support of the intermediate scattering 
function by including a step function giving 

Fs{k, <) = [1 - (1 - q) DkH] e (1 - (1 - 9) DkH) , 

(16) 

where Q {y) = 1 for y > 1 and O (y) = otherwise. 
Substitution of H16|l into Ea. H12|l shows that (|16|l satisfies 
the (j-diffusion equation since the time-derivative acting 
on the step function gives a Dirac delta function with 
a prefactor that goes to zero. Note that in this case, 
since the intermediate scattering function is a function 
of Dk^t, it immediately follows that in real space it is a 
function of r^/Dt so that (r^) ~ Dt which corresponds 
to normal diffusion, albeit with non-classical distribution 
functions. 

In one dimension, inversion of the Fourier transform 
(115(1 in the case q > I gives 



G,>i{r,t) 



9-1 



^^iq-l)Dt 



K 



(17) 

where x = , ^ , and K^, (x) is the modified Bessel 

function of the second kind. The structure of this result 
should not be too surprising since we know that the gen- 
eral solution to the classical diffusion problem is given in 
terms of Bessel functions jS] . For some rational values of 
q it is possible to express the Bessel function in terms of 
elementary functions in which case the solution assumes 
the form p{x) exp(— x) where p{x) is a polynomial in x. 
For instance for q = 3/2, we obtain 



G3/2(r,i) 



1 



■ {x + l) exp (— .t) 



(18) 



Figure 1 shows a numerical evaluation of the solutions 
for various values of g, and, as q increases, it does indeed 
appear to decay as an exponential. 

We now show that H17|l does give the correct Gaussian 
limit as q 1. We make explicit the dependence on q 
by writing 



i(a;) = ^ , , 

5-1 2 \^q-\ 



(19) 



where r* = r/V Dt. Introducing ly = — i in H17|l 
gives 



Gg>i(r,t) 



X K, 




(20) 
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Using the asymptotic expansion for v ^ 1 



1 



X exp (^-1^ + , 
with z = r* V^^iZ^ gives, for fixed r*, 



X exp ( J/ In 2 + -jf*'^ I exp ( — —r 



2v 



wherefrom we obtain 



1 *2 

exp I —V — —r 



/2Dt 



V exp 



and, using Stiriing's formula, we find 



G,>i(r,i) 



1 



— . exp I I , 



(21) 



(22) 



as expected. Fig. 2 clearly shows that the distribution 
approaches a Gaussian in the limit v \ ox q —* \. 
For q < 1, Fourier inversion of (|f 5|l yields 



G,<i(r,t) 



with X = 



1-9 



1 ...-^ 



20F^(f-g) Bt 



(23) 

, , and where J^, is the Bessel function 

of the first kind. Here also, for some values of q, this 
result can be written in terms of elementary functions in 
the form l(x) cosx + m{x) sinx where l(x) and m{x) are 
polynomials in x. As an example, for q = 1/2, we have 



Gi/2(r,i) = 



TT 



-1 



X 

2Dt V2 



(3 sin: 



3a; cos a; — x sin a;) 



(24) 

Figure 2 shows that as g ^ 1 the distribution approaches 
a Gaussian, and Fig. 3 shows Gq{r,t) for several values 
of g < 1. In this case, we observe that the distribu- 
tion goes sometimes negative (for small values of q) and, 
hence, becomes unphysical. Note that after the distri- 
bution first becomes negative, its magnitude is always 
very small so that in some sense it appears to be ap- 
proximating a function with finite support. Therefore, 
we must either (a) restrict the range of allowable val- 
ues of q to g > 1, or (b) allow values of g < 1 with 
the understanding that the distribution is an approxima- 
tion to a positive-definite function, perhaps one having 



finite support. It is worth noting in passing that option 
(b) is quite commonly used in other areas of statistical 
mechancs. For example, the Chapman-Enskog solution 
to the Boltzmann equation results, at all orders, in a 
distribution which is not positive-definite [lo| . 

It is interesting to ask what conventional diffusion pro- 
cess would yield the above results. Casting the general 
solutions H17(l and (|23(l into the form of a generalized 
diffusion equation in real space 



dt 



Dgir,t) 



gives, for g > 1, 



i?,>i(r,i) = D 



g-1 N{r,t) 
2 M{r,t) 



(25) 



where, with the notation v = — i, 
N{r,t) 



q- 



1 - 2j/, 



M{r,t) = 1 
and, for g < 1 
D. 



(i-ll^(l-2.)(iV(M) + l) 



[r,t) = D 



1-g N'{r,t) 
2 M'{r,t) 



(26) 



where, with the notation v 



N'{r,t) 



l-q ^ 2' 



M'{r,t) = -l + il-|l^(l + 2zy)(iV'(r,i)-f 1) . 

Two typical examples (for q = 3/2 and q — 1/2) of the 
generalized diffusion coefficient are |ll| 



^3/2 (r,t) = — 



where D' = ^D, and 



1 



(27) 



/D't . 



with 



D' x^ sin a; (3 — a;)(xcos x — sin x) 



4 sin X - 



l)(cos X 



(28) 



ID't 



Although the diffusion coefficient Dq exhibits a gener- 
alized form because of its dependence on the distribution 
function Gq (r, t) , diffusion remains normal in the sense 
that -^^^ scaling is preserved as we already observed and 



4 



as can be seen even more explicitly by computing the 
mean square displacement 

/■oo 

(rHt)) = / drr'G,{r,t), (29) 
Jo 

with Gq{r,t) given either by H17|l or (|23|l . which, with 
r* — r/\^Dt, yields 

/•oo 

{r^{t)) = Dt2 / dr*r*^Gq{r*) = 2Dt. (30) 

In this paper, we have shown how a standard deriva- 
tion of the diffusion equation can be generalized into a 
" nonextensive" form without the introduction of ad hoc 
modifications. Instead, we simply generalize the ansatz 
used to resum the expansion of the intermediate scat- 
tering function. For an intial condition of a Dirac delta 
function in real space, the solutions of the resulting dif- 
fusion equation in Fourier space take the form of the q- 
exponential which commonly arises in nonextensive sta- 
tistical mechanics. This solution corresponds to normal 
diffusion, but with non-classical distributions. There- 
fore the g-generating function approach developed here 
should be useful to similarly generalize other derivations 
of diffusive transport. 



FIGURE CAPTIONS 

FIGURE 1 : Distribution Gq>i(r, <) as a function of 
r* = r/Voi. 

FIGURE 2 : Distribution Gq (r, t) as a function of 
r* = r/^/Dt for q — 1.05 and q = 0.95 compared to 
the Gaussian solution. 

FIGURE 3 : Distribution Gg<i(r, i) as a function of 
r* = r/y/m. 
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